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Abstract 

We discuss the scale dependence of a skewed parton distribution of the pion ob- 
tained from a generahzed hght-cone wave function overlap formula. Using a simple 
ansatz for the transverse momentum dependence of the light-cone wave function 
and restricting ourselves to the case of a zero skewedness parameter the skewed 
parton distribution can be expressed through an ordinary parton distribution 
multiplied by an exponential function. Matching the generalized and ordinary 
DGLAP evolution equations of the skewed and ordinary parton distributions, re- 
spectively, we derive a constraint for the scale dependence of the transverse size 
parameter which describes the width of the pion wave function in transverse mo- 
mentum space. This constraint has implications for the Fock state probability 
and valence quark distribution. We apply our results to the pion form factor. 



Skewed parton distributions (SPDs) provide a link between exclusive and inclusive 
quantities of QCD 0, ^. Among some of their well known properties is their relation to 
ordinary parton distributions and hadronic form factors via so-called reduction formulas. 
Moreover, their evolution behaviour has been investigated and generalized DGLAP evo- 
lution equations have been derived. Only few is known, however, about their particular 
form and so for applications to physical processes one has to resort to specific models. The 
authors of discussed SPDs in the context of soft contributions to large angle Compton 
scattering and form factors and proposed a generalized Drell-Yan overlap formula |^ for 
SPDs in terms of light-cone wave functions (LCWFs). It was shown that in a special kine- 
matical region, where the plus component of the momentum transfer and consequently the 
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skewedness parameter vanishes, and with a Gaussian ansatz for the transverse momen- 
tum dependence of the LCWFs, SPDs can be expressed by ordinary parton distribution 
functions multiphed by an exponential t dependence. 

In the present paper we will investigate the consequences of the evolution equations for 
this phenomenological model of SPDs. As we will argue below, the combined evolution 
of the SPD and the ordinary parton distribution enforces a condition upon the scale de- 
pendence of the transverse size parameter. This parameter appears in the SPD through 
the transverse momentum dependence of the LCWFs and describes the width of the wave 
function in /c^-space. We will derive a model dependent intregro-differential equation for 
the scale dependence of the transverse size parameter, which we will solve numerically. 
Moreover, we will show that the scale dependent transverse size parameter induces a scale 
dependence of the Fock state probability of the lowest Fock state and the corresponding 
valence distribution. As an application to a physical process it is natural to consider soft 
overlap contributions to the pion form factor. We conclude with our summary. 

In this work we use the conventions of Radyushkin and denote the SPD of a parton 
with flavour a in the pion by It is defined by a bilocal matrix element of quark 

field operators: 

where the notation tpa{z~) indicates that the argument of the operator has vanishing light- 
cone plus and transverse components. The following reduction formulas [0, Q are general 
properties of SPDs. In the forward case one regains ordinary parton distributions: 

^^^,{x,t = 0) = qaix), (2) 

and by integrating the valence distribution one obtains the pion form factor: 

F^{t) = dx^^{x,t), (3) 

where we have defined ^^{x, t) = ^^{x, t) — t) and used e„ — = 1. 

Our starting point is the overlap formula for SPDs in the region ( < x < 1 0]: 

t) = ^ 1 E / [dx]N[d'k,_U S{x - xi)^%^{x'„ k^J ^^^{x., k^O, (4) 

-'- ^/^ 1,13 

where denotes a particular Fock state, the summation index / runs over all quarks of 
flavour a and f3 refers to different spin flavour combinations of partons in a given A^-particle 
Fock state. The integration measures are defined by 



TV . . 

[dx]N = Y[dXnSil-Y,Xmj, 

n=l ^ ni ^ 



[d'k.u - (Y^;^^ n ^^k,„5(^)(x:k.„) (5) 
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and the arguments of the initial and final wave function are related by 



X 



I 



k'^,. = k^, + — -^A^, (6) 



with j being the index of the active parton and i being the index of the spectator partons. 
A_L is the transverse momentum transfer between the initial and final hadron. 

Following the authors of ||, 0], we write the soft A^-particle LCWF of the pion 
in terms of the distribution amplitude (pj^p and a transverse momentum dependent part, 
for which we make a Gaussian ansatz: 



^Ni3{xi, k_Li) = MNp(t>Ni3{xi) ^^^'^ exp - ^ — 



2^2 \N-1 



X1X2 .. .Xn 



^ k^ 



=1 ■^i 



(7) 



where J^n/b is a normalization constant and, for obvious reasons, ajv is called the transverse 
size parameter of the iV-particle Fock state. 

As we have already mentioned we will restrict ourselves to the special case of a zero 
skewedness parameter, i.e. ( = 0. It has been shown in [Q that the Gaussian k^- 
dependence then leads to the following simple representation of an SPD in terms of an 
ordinary parton distribution function multiplied by an exponential: 



^fj',\x,t)=qi''\x) exp 



2 X 



where we use the common notation x = 1 — x. The origin of the appearance of the parton 
distribution in (Rf) is its relation to LCWFs via the expression [H] 



J2 [[dx]N[d'^^±]NSix- Xl)\^Nf3iXi,'k±^)\'^. (9) 
1,(3 •' 



As will become obvious immediately it is useful to approximate the transverse size param- 
eters of all Fock states by a common value a^r, i.e. we set 

~ Oat for all A^. (10) 

This is in general a rather rough approximation. However, from the exponential in Eq. (|) 
together with the fact that with increasing A^ the functions qa^\x) are proportional to 
increasing powers of (1 — x) it is clear that large x dominate at large t and only a few 
of the lowest Fock states contribute to phenomenological applications. We can now sum 
over all Fock states, 

•^C=o(^,t)=E-^cloHa:,t), (11) 

N 

SO as to obtain a representation in terms of the full quark distribution function. Since 
we will discuss the pion form factor later on and in order to avoid complications in our 
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discussion of evolution by quark-gluon mixing we consider the valence distribution from 
now on: 



^l=Q{x,t) = q^,{x) exp 



1 9 x ■ 



(12) 



The SPD (0) is completely independent of the particular form of the pion distribution 
amplitude. Thus, we need not to specify the pion distribution amplitudes 0iv/3 in A 
corresponding expression for the SPD of the nucleon has been suggested in Refs. [§, |] and 
in case of the pion in Ref. |T^ and also recently by the authors of |1T 



At this stage, the SPD (|T2D depends on a scale iJ? only through the ordinary parton 
distribution. However, as we are going to show the evolution equation of the SPD forces 
the transverse size parameter a^r to be scale dependent as well. The scale dependence of the 
l.h.s. of expression (|T2|) is described in terms of a generalized DGLAP evolution equation, 
where for C = the modified evolution kernels are known to reduce to the ordinary DGLAP 
kernels 

/.^A^^o(-,t;;.^) = ^ (13) 

The r.h.s. of (|12D is given in terms of a parton distribution obeying ordinary DGLAP 
evolution, multiplied by an exponential. Obviously, both the evolution equations of the 
SPD and of the parton distribution can only be fulfilled simultaneously if 



(14) 



so that for the r.h.s. of ([T^) we have: 



d 
exp 



qv{x;^^) exp 



2tt 



dy 

X y 



V?/ / AT. 



2x 



(15) 



Equating expressions (|I3|) and (|T5D and resolving for the derivative of a7r(/^^) we find an 
integro-differential equation for the scale dependence of the transverse size parameter: 



/i 



~~d^ 



1 X 
2x 



1 -1 



tqv{x; /i 



X 



as(/i^ 
27r 



/'*P„(£)(l_exp[iaV)(^-2)< 

Jx y ^y^ \ 12 ^y x^ . 



(16) 



This equation is a consequence of our particular model ([l2|) . As we can see immediately, 
the evolution of the transverse size parameter is driven by the difference of the SPD and 
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the ordinary parton distribution. This equation may be solved numerically using an iter- 
ative method similar to the one employed for the numerical solution of ordinary DGLAP 



equations, see for instance |12]. In order to determine the initial value of we consider 



the LCWF of the lowest Fock state. It is commonly accepted that the form of the pion's 



two-particle distribution amplitude is close to the asymptotic one |T3[, (l)^{x) = Qxx, to 
which we will restrict ourselves in the following. The parameters of the two-particle LCWF 
are then completely fixed from various decay processes p|. The normalization follows from 
the 7r+ /^^t'/i decay and it is given by J\f2 = /^/(2a/6), where ~ 131 MeV is the well 
known pion decay constant. The two-photon decay of the uncharged pion, 7r° — 77, leads 
to a constraint for the transverse size parameter, for which a value of = 0.86 GeV~^ 
is found. As the corresponding scale we choose /Xq = 0.25 GeV^ which is a typical scale 
for light mesons [1^. From now on we will employ the GRS parametrization of the 
valence quark distribution 

It is interesting to note that the singularity of the splitting function is canceled by the 
expression in curly brackets of Eq. ([I6|) . As we can see further, there are two free variables 
in Eq. (|TBp: the longitudinal momentum fraction x and the momentum transfer t. This 
means that our approach does not a priori exclude an explicit x- and t-dependence of a^^. 
Intuitively, however, we expect transverse quantities such as the mean square transverse 
momentum (k^^) neither to depend on the longitudinal momentum nor on the momentum 
transfer. The numerical investigation indeed only shows a very mild variation of a7r(/i^) at 
different values of large x and t. If this were not the case the ansatz (|1^) would have to be 
abandoned. For definiteness, in what follows we will choose x in Eq. ([T6| ) to be the mean 
value of the momentum fraction contributing to the pion form factor which is given by 

1 /•! 

^^)* •= 7^ / (17) 
K{t) Jo 

Taking —t = 10 GeV^ we find a value of {x)t = 0.75. Varying —t between 1 GeV^ and 
10 GeV^ a.7r(yU^) only changes by less than 4%. 

The scale dependence of the transverse size parameter is shown in Fig. |l[ As can be 
expected a-j^ depends only moderately on /i^ and the decrease is weakened with increasing 
scale. In Fig. |^ we have plotted the valence SPD (|12D at two different scales. In order 
to see which quantitative effect the scale dependence of has on the SPD, we compare 
J^l=o{.x,t] n'^) with frozen (dashed line) and running a.„ (solid line), respectively, at jj^ = 
10 GeV^, where the scale dependence of the transverse size parameter causes a shift of the 
SPD towards smaller momentum fractions. At = 0.25 GeV^ both curves coincide by 
definition (dot-dashed line). 

At constant a,r the parton distribution of the lowest Fock state which results from the 
LCWF (0) is scale independent. With increasing scale, however, one expects a damping of 
the parton distributions of lower Fock states since an increasing number of virtual quark- 
antiquark pairs and gluons is produced so that higher and higher Fock states become 
occupied. Quantitatively, this damping effect emerges through the Fock state probability 



5 



Scale dependence of a„ 
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Figure 1: Scale dependence of the transverse size parameter a^r with an initial value of 
= 0.86 GeV~^ at fi^ = 0.25 GeV^. As described in the text we have chosen x = {x)t = 
0.75. The two curves demonstrate the weak dependence on t. 



in our approach. The Fock state probability of the A^-th Fock state is defined by 



The running a^^ induces a scale dependent Fock state probability. As we have already 
discussed below Eq. (p!6| ) the normalization of the LCWF of the pion's lowest Fock state 
is fixed so that for P2 we have explicitly 



which coincides with the well known value of P2 = 0.25 at fiQ = 0.25 GeV^ due to our 
choice of the initial value of a^^. As we can see immediately the Fock state probability of 
the lowest Fock state decreases with increasing scale since it is proportional to a^(yU^). The 
value of P2 at fi"^ = 100 GeV^ is reduced to about 0.14. 

We can now write down the corresponding scale dependent valence quark distribution. 
Using the asymptotic form of the pion's two-particle distribution amplitude in expres- 
sion we obtain 

u^^\x; fx^) = 6 P2ifi'') XX, (20) 

which is shown in Fig. |^ at two different scales, where we also compare with the GRS 
parametrization. The plot clearly shows the anticipated damping of the valence distribution 
with increasing scale. As discussed below Eq. (^) only a few Fock states contribute to 
the parton distribution at large x. We thus expect the valence distribution u^\x) of 
the lowest Fock state to approximate well the full valence distribution Uy{x) at large x. 
With a fixed transverse size parameter, Fig. |^ shows that this expectation is not fulfilled 
at large scales for x larger than 0.6 since Uy{x) is shifted to smaller values of x while 




(18) 



P2i^^')=2n^f:al{^^% 



(19) 
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Valence SPD 
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Figure 2: Evolution of the valence SPD (0) with constant and running a^. We use the 
GRS parametrization of the forward valence distribution as stated in the text. At the 
starting scale /ig = 0.25 GeV^ both curves naturally coincide. The momentum transfer —t 
is arbitrarily set to 1 GeV^. 



u^y\x) remains constant. We see that switching on the scale dependence of a^, complies 
with our expectation provided that we take into account the theoretical and experimental 
uncertainties of the analysis of the Drell-Yan process 7r~p — > fi~^fi~X, from which the GRS 
parametrization of the valence distribution is extracted. In particular, for a; > 0.75 the 
GRS valence distribution is less reliable since in that region the parametrization of the 
proton structure function, which is used as an input in the analysis, is an extrapolation. 
Our approach provides a clear improvement compared to the results of Ref. ||16|, where a 



constant transverse size parameter has been used and where thus a comparison of the full 
and the N = 2 parton distributions has been possible only at low scale. 

We will now consider the pion form factor, given by Eq. @, as an application.^ Since 
the relevant scale in the overlap formula is given by the momentum transfer t = — with 

being the photon virtuality in a frame where the plus component of the momentum 
transfer vanishes, such that C = 0, it is natural to identify /i^ with —t. The result is shown 
in Fig. ^ (solid line). For comparison we also plot the form factor neglecting the evolution 
of both and a,r (dashed line), and with fixed a^, and running (dot-dashed line), 
respectively. The three curves show that the scale dependence of a-,, roughly compensates 
the evolution of the GRS valence distribution in the region 5 GeV^ < — t < 10 GeV^. 
For —t<5 GeV^ the running a^r even provides a slight enhancement of the theoretical 
prediction. 



The authors of H, IT6| considered the contributions of the lowest Fock state only, without 



^As can easily be shown numerically the ansatz (^ of the pion's LCWF does not provide significant 
contributions to the Drell-Yan overlap formula of the pion form factor for a; 1 and fcj^ — )■ and so 
contrary to what has been claimed in for instance, the LCWF (^ does represent well soft QCD 
contributions. 
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Valence distribution 
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Figure 3: The scale dependent valence quark distribution pO| ) of the lowest (A^ = 2) Fock 
state (long dashed and solid lines) compared with the GRS parametrization (dashed and 
dot-dashed lines). As we discuss in the text to the right of the bar this parametrization 
becomes less reliable. 



taking into account evolution effects of the transverse size parameter. Comparing their 
results with the dashed line of Fig. ^, where we have frozen the evolution scale of both the 
transverse size parameter and the valence distribution, we see that their prediction stays 
below ours. This has to be expected since we take into account the contributions of all 
Fock states. 



The factorized GRS ansatz of Ref. |TT[ is essentially identical to our expression ([1^ ) 



again with a scale independent transverse size parameter and the GRS valence distribution 
at a fixed scale. Their corresponding prediction of the pion form factor is somewhat higher 
than ours. This is due to the fact that the parameter Aq used in [|TT| corresponds to a 
transverse size parameter of a^^ ~ 0.77 GeV~^, which is smaller than the value used in the 
present work. 

The predictions of the hard contributions to — tF^(t) alone, ranging from 0.08 GeV^ in 
T6| to about 0.16 GeV^ in |jl9|, obviously cannot account for the experimental data. We 



would like to point out that the sum of the hard part of Ref. [0 and our prediction of the 
soft part is in very good agreement with the new data between 0.6 and 1.6 GeV^ presented 
in Ref. [^. Note that in Ref. strong cancellations between soft parts and hard parts 
of higher twist have been found, leaving small non-perturbative contributions. 

To summarize, we have shown that the evolution equations for SPDs leads to a further 
constraint for phenomenological models of SPDs which are expressed through ordinary 
parton distributions. Starting from the generalized Drell-Yan formula, where the SPD for 
the A^-th Fock state of the pion is written in terms of an overlap integral of A^-particle 
light-cone wave functions, and specializing to the case of a zero skewedness parameter the 
SPD equals a Fock state parton distribution multiplied by an exponential function. Making 
further simplifications by assuming a common transverse size parameter for all Fock states 
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Pion form factor 
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Figure 4: The pion form factor with constant and running a^^. The dashed hue shows the 
case of constant a^r = 0.86 GeV~^, where we have also fixed the scale of the GRS valence 
distribution at 1 GeV^. Data are taken from Ref. IIT? . 



we have obtained the full SPD, to which we have then applied the evolution equations. 
We have matched the generalized and ordinary DGLAP equations for the skewed and 
forward parton distributions, respectively, which has resulted in a constraint for the scale 
dependence of the transverse size parameter. This in turn has led to a scale dependent 
Fock state probability and valence quark distribution of the lowest Fock state of the pion. 
The application to the soft overlap contributions to the pion form factor has shown a slight 
enhancement of the theoretical prediction in the few GeV^ region, which is in complete 
agreement with new data. 

Finally, we would like to remark that LCWFs of the form (^, which are modified by an 
effective mass term, have also been discussed in the literature, see, for instance IHl and first 



Ref. of [1^. However, our discussion of the evolution effects shows that in principle one 
has to take into account the scale dependence of the effective mass as well, which would 
reduce the influence of a mass term in the LCWF with increasing scale. 
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